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We present a study of "nodal semimetal" phases, in which non-degenerate conduction and va- 
lence bands touch at points (the "Weyl semimetal") or lines (the "line node semimetal") in three- 
dimensional momentum space. We discuss a general approach to such states by perturbation of the 
critical point between a normal insulator (NI) and a topological insulator (TI), breaking either time 
reversal (TR) or inversion symmetry. We give an explicit model realization of both types of states 
in a NI-TI superlattice structure with broken TR symmetry. Both the Weyl and the line-node 
semimetals are characterized by topologically-protected surface states, although in the line-node 
case some additional symmetries must be imposed to retain this topological protection. The edge 
states have the form of "Fermi arcs" in the case of the Weyl semimetal: these are chiral gapless 
edge states, which exist in a finite region in momentum space, determined by the momentum-space 
separation of the bulk Weyl nodes. The chiral character of the edge states leads to a finite Hall 
conductivity. In contrast, the edge states of the line-node semimetal are "fiat bands": these states 
are approximately dispersionless in a subset of the two-dimensional edge Brillouin zone, given by 
the projection of the line node onto the plane of the edge. We discuss unusual transport properties 
of the nodal semimetals, and in particular point out quantum critical-like scaling of the DC and 
optical conductivity of the Weyl semimetal, and similarities to the conductivity of graphene in the 
line node case. 



I. INTRODUCTION 



The study of systems, distinguished by topology rather 
than symmetry, is an increasingly important theme in 
modern condensed matter physics. This paradigm shift 
has gained more momentum recently, with the discov- 
ery of the time-reversal (TR) invariant topological insu- 
lator (TI).^ ^5 Apart from reinvigorating the interest in 
topological phenomena in solids generally, this discovery 
has drawn particular attention to the momentum-space 
topology of the electronic band structure of solid crys- 
talline materials, or, in a more general context, to the 
momentum-space topology of fermionic ground states. 
The most common view of topologically nontrivial elec- 
tronic phases is that these are states of matter, which are 
insulators in the bulk, yet have metallic edge or surface 
states, which result from the nontrivial momentum-space 
topology of the bulk band structure. The gap in the 
bulk electronic spectrum is what makes the topological 
ground state insensitive to small perturbations and pro- 
tects (perhaps in combination with a discrete symmetry, 
such as TR) the metallic surface states. The appearance 
of such a robust metallic surface state in a bulk insula- 
tor is the main experimentally-observable manifestation 
of topological order. This, however, is an oversimplified 
view. Very recent work has shown that certain special 
types of gapless band structures can in fact also be topo- 
logically nontrivial and give rise to robust gapless surface 
states.lS"^ Gapless topologically-nontrivial band struc- 
tures are characterized by the presence of point or line 
nodes, i.e. points or lines in the three-dimensional (3D) 
momentum space, at which two distinct bands touch each 
other accidentally. While such accidental band-touchings 
have been known to exist and studied since the early days 



of the band theory of solidSj^^^ only much more recently 
have their nontrivial topological properties been noticed 
and their significance appreciated, st artin g, in particular, 
with the pioneering work of Volovik.lSE^ 

Topological properties of the accidental band-touching 
nodes depend crucially on their co-dimension.'^^ The 
point nodes, which have an odd co-dimension 3 = 3 — 
0, are the most robust variety. The band structure 
near such point nodes is described by a massless two- 
component Dirac (Weyl) Hamiltonian, and is topolog- 
ically equivalent to a hedgehog in momentum space.131 
The only way to eliminate such a momentum-space 
hedgehog is to annihilate it with an antihedgehog, i.e. 
Dirac point of opposite chirality. The theorem of Nielsen 
and NinomiyeP^ guarantees that Weyl nodes always oc- 
cur in pairs of opposite chirality. Such pairs are topolog- 
ically stable if the opposite-chirality partners are sepa- 
rated in momentum space, thus precluding their mutual 
annihilation. When both TR and inversion symmetry are 
present, however, energy bands are two- fold degenerate 
at all momenta, and non-trivial contact is between pairs 
of bands. This has vanishing probability at generic points 
in momentum space, and occurs only by tuning one pa- 
rameter at special time-reversal invariant momenta. In 
the latter case, the four bands crossing can be viewed as 
a pair of opposite chirality Weyl nodes which occur at the 
same point in momentum space, and are thus not stable 
- hence the need for a tuning parameter. The separation 
of the two Wey l no des in momentum space is achieved 
by breaking TFiP^i or inversion symmetry. The result- 
ing topological Weyl semim etal p hase can then be shown 
to possess chiral edge states^I^ (this is the only known 
example of a state with topological chiral edge states 
which is intrinsically three dimensional) and a nonzero 
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Hall conductivity, proportional t o the separation of the 
Dirac nodes in momentum spaceP^I^ 

A line node has an even co-dimension 2 = 3 — 1 
and does not possess the absolute topological stability 
of a point node (this even-odd dichotomy is an exam- 
ple of Bott periodicity and can be understood within 
-ftT-theor jEl) . However, imposing certain discrete sym- 
metries can stabilize line nodes, i.e. they may be stable 
with respect to all pe rturb ations, not violating a spe- 
cific discrete symmetry.f^^SI As for point band touching, 
imposition of both TR and inversion symmetries is too 
restrictive for line touchings to occur. Thus we consider 
here the case where TR is broken. However, the role of 
symmetry in the case of the line nodes is more complex 
than for Weyl nodes, and we will see that, while the line 
contact between bands can be stabilized, this requires a 
delicate, but physically achievable, combination of dis- 
crete symmetries (other than TR). Moreover, the condi- 
tions, which stabilize a line contact between two bands, 
are distinct from those which force this line contact to 
have constant energy (and are in general insufficient to 
guarantee constant energy). Thus we will argue that by 
physical symmetries a line node cannot be stabilized at 
the Fermi level. Interestingly, even a line contact, which 
is not at constant energy (which we will continue to call 
a line node though it is a slight abuse of terminology in 
this case), has topological properties and can be related 
to the surface spectrum. The distinguishing character- 
istic of these surface states is that they exist inside the 
"direct gap" between conduction and valence bands in a 
finite area in the 2D surface Brillouin zone (BZ), whose 
boundary is determined by the projection of the bulk 
nodal line onto the plane of the surface. In the sim- 
plest models, which may be a reasonable approximation 
in some cases, these surface states are entirely dispersion- 
less, i.e. form a "topological" flat band.^^ 

In this paper we discuss a particularly simple realiza- 
tion of both a Weyl and a line-node semimetal phases, 
which obtains in a multilayer heterostructure material, 
composed of alternating layers of a TI and a normal insu- 
lator (NI) material.'Sl This can be understood as a simple 
means of constructing an underlying 4-component Dirac 
point with both TR and inversion symmetry, describing 
the (unstable) critical point between a bulk NI and TI 
phase. Then by a judicious choice of TR breaking per- 
turbation, we can realize both stable nodal phases (the 
alternative case, in which inversion symmetry is broken, 
has recently been discussed in Ref . [T2l) . Here we imagine 
doping with magnetic impurities, which are presumed to 
ferromagnetically align, but applying an external mag- 
netic field will also serve, though this introduces some 
changes in low energy properties due to the influence of 
the orbital component of the fleld. We provide a char- 
acterization of the edge states in the two topological 
semimetal phases and also discuss their transport and 
optical properties, which are unusual and should be a 
focus for experimental studies. Our results for the trans- 
port and optical properties of the Weyl semimetal apply 



equally well to proposed bulk realizations of the state, 
such as in pyrochlore iridates.® 

The paper is organized as follows. In Sec. |TT] we dis- 
cuss the general theory of perturbed 4-component Dirac 
points, and classify the perturbations, which give rise to 
point nodes, line nodes, and Fermi surface phases. These 
results may apply very generally, and allow design of any 
desired nodal state, once a physical meaning, appropri- 
ate to a speciflc system, is given to each of the Dirac 
matrices. In section |III| we review the point-node, or 
Weyl, semimetal, discussed previously in Ref. [HI in this 
context, and discuss the effects of an orbital magnetic 



fleld on this state. In section IV we give a detailed dis- 
cussion of the line-node semimetal, including its unusual 
"flat-band" surface states and effects of an orbital fleld. 
Section [V] discusses the conductivity of both the Weyl 
and line node states. We conclude in section IVII with a 
brief summary and discussion of our results. 



II. GENERAL THEORY: PERTURBED DIRAC 
POINTS 

In this section, we consider generally how semimetallic 
phases may emerge from perturbations of a system close 
to a TI to NI transition in a nearly time-reversal and in- 
version symmetric system. The result is that point node, 
line node, and metallic Fermi surface states are possible, 
depending upon the nature of the perturbation. 



A. Dirac equation and matrices 

When both TR and inversion (I) symmetry are present, 
a direct NI-TI transition is possible by the formation of a 
massless 3-1-1 dimensional four component Dirac fermion 
at the critical point .'^ This occurs at a time- reversal in- 
variant momentum, which for simplicity we take to be at 
the F point k = 0. The k ■ p expansion about the F point 
then generically takes the form: 



Ho 



(1) 



where 7^ (/i = 1 • • • 5) are the five 4x4 Dirac matrices 
in an appropriate basis. We have rescaled coordinates to 
set the Dirac velocity to unity for simplicity. Since the 
momentum components ka are odd under both TR and 
I, 7i,72, and 73 must also be odd under both TR and I, 
while 74 is even under both. Then, since 75 may be ob- 
tained as the product of the four other gamma matrices, 
it is odd under both TR and I. The full space of Her- 
mitian 4x4 Hamiltonians is spanned by including the 
identity and another 10 matrices, 706 = — |[7a,7fc] with 
a < b. One can deduce their transformation properties 
from those of the 7^. The gamma matrices can be sep- 
arated into three vectors, b,b',p and one scalar e, with 
transformation properties given in Table [ij When both 
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operator 


TR 


I 


b 


= (723, 7i3, 712) 


-1 


+1 


P 


= (714,724,734) 


+1 


-1 


b' 


= (715,725,735) 


-1 


+1 




£ = 745 


+1 


-1 



TABLE I. Transformation properties of Dirac operators. 



TR and I symmetry are preserved, all 10 of these ma- 
trices are prohibited from entering the Hamiltonian by 
symmetry. Only the mass m is allowed, and there is thus 
a single tuning parameter to access the massless Dirac 
point, which separates the TI and NI phases. 



3. u ■ V = perturbation 

When both u and v arc non-zero, the spectrum de- 
pends upon their relative angle. When u • v = 0, 
it can still be diagonalized analytically. Taking u = 
(«,0,0),v= (0,z;,0), 



ei(u-v==0) ==± 



1/2 



(5) 



From Eq. ([5|, one finds that when u^+v^ > m^, there are 
two Weyl nodes at = = 0, fci = ±\/v? + v"^ — m?, 
and when v? +v'^ < m^, there is a full gap. 



B. TR breaking perturbations 



4- u II V perturbation 



Let us now consider what happens to this critical point 
when either TR or I symmetry is relaxed. First consider 
relaxing TR, but preserving L In this case, the terms b 
and b' may be added, and the most general Hamiltonian, 
which involves constant coefficients perturbing the Dirac 
point is: 



•b'. 



(2) 



In general, this is too difficult to diagonalize analytically. 
For several simple cases, however, it is possible. 



1. b perturbation 

For V = 0, by an 0(3) rotation, we may choose the 
Hamiltonian in the form Hi = Hq + uhi . This gives the 
spectrum: 



ei(«-0)-±^|^ 



2 + fc2±u +kl + kl. (3) 



This gives two stable Weyl nodes with ei = when |u| > 
m, with fci — zbVw^ — TTt^ and k2 — k^ — 0. 



For u • V 7^ 0, in general Hi cannot be diagonalized 
analytically. An exception is the case to = and u || v, in 
which case we can take, e.g. u = (m, 0, 0) and v = (w, 0, 0) 
by an 0(3) rotation, and the spectrum is: 



ei(TO = 0) = ± 



+v'^ + k^ 



±2 



From Eq. ([6]), if |u| > |t>|, one has two Weyl points at 
fci — ±^ — v-^, ^2 = ^3 = 0, while for |u| < \v\, there 
is a ring node at fci =0, fc| -I- fc| — v'^ — u^. While we can 
no longer find the spectrum analytically when to 7^ 0, 
in this case, we find numerically that the response to 
such a mass is distinct from the situations above. For 
small TO 7^ 0, the point and line nodes expand into Fermi 
surfaces: small pockets for |u| > l^l, and a torus for 
|u| < |z;|. As TO is increased, these surfaces evolve and 
eventually shrink to a point at some threshold to* , above 
which there is again a gap. Thus in this case neither 
point nor ring nodal states are stable, and instead the 
NI-TI transition is converted to an intermediate metallic 
state. 



C. I breaking perturbations 



2. b' perturbation 

Next consider the case u = 0. In this case, by a similar 
rotation, we have Hi = Hq + vb'i, and: 




2+fc2+fc2±„ +fc2. (4) 



ei{u = 0) = ± 



In this case, when |i;| > to, there are two bands whose 
energies touch along a circle, defined by fc|-|-fc| = v^—m?, 
fci = 0. 



If inversion symmetry is broken, but time reversal is 
preserved, the most general Hamiltonian with constant 
coefficients is of the form: 



H2 ~ Hq + ■ ^ + \e. 



(7) 



Without loss of generality, we can use an 0(3) rotation 
to choose w = (ii;,0,0), and the resulting Hamiltonian 
can be diagonalized to obtain: 



£2 =± 



rri^ + w'^ + \^ + k^ 



±2y^A2fc2 + (u;2 + A2)(fc 



1/2 



(8) 
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Here the spectrum is fully gapped whenever m ^ Q. For 
m = 0, there is a gapless nodal line located at ki = 0, 
/cg — w'^ + \^ . Thus at this level of approximation, there 
remains a direct NI-TI transition when m = 0, but with 
a critical nodal line formed at the transition. In fact, this 
is an artifact of the approximation we have made, that 
the coefficients to. A, w are momentum-independent. As 
shown in Ref. 12, when proper momentum dependence 
is included, this transition point broadens into a Weyl 
semimetal phase, with a minimum of 4 nodal points. 



III. POINT-NODE (WEYL) SEMIMETAL IN A 
TI MULTILAYER 

A. Model and connection to Dirac equation 

In the previous section, we observed that a point node 
state could be generated by certain time-reversal sym- 
metry breaking perturbations of the TI-NI Dirac critical 
point (e.g. Sec. II B 1). Here we discuss the specific case 
of this mechanism in a model of a TI multilayer het- 
erostructure, introduced by two of us in Ref. IHl 



Here z, j label individual TI layers, separated by NI spac- 
ers, Ag is the tunneling matrix element between the top 
and bottom surfaces within the same TI layer, A^ is the 
tunneling matrix element between the top and bottom 
surfaces of nearest-neighbor TI layers, and kj^ is the mo- 
mentum in the 2D surface BZ of each TI layer. Without 
loss of generality we will assume that A5, A^i > 0. Such 
a multilayer structure exhibits a critical point between a 
strong 3D TI, when A^ > Ag and an ordinary 3D insu- 
lator when A5 > A^i. The critical point, A5 — A^i, 
realizes the 4-component Dirac fermion, which is the 
starting point of the previous section. It occurs here 
when the gap vanishes at a single point in the 3D BZ 
kx = ky = 0, kz = Tr/d, where d is the superlattice period 
of the multilayer. The momentum-space Hamiltonian, 
expanded to leading nontrivial order near this point, is 
given by: 



(10) 



where vp = d^/AsAj:). This is the Hamiltonian of a 4- 
component massless Dirac fermion, equivalent to Eq. ([T]) 
in Sec. |ll] after a rescaling of coordinates k^ — >■ k^/vp, 
ky — >■ ky/vp, kz — >■ kz/vp- We can identify from it a 
specific physical realization of the first three gamma ma- 
trices: 



7i 



72 



73 



(11) 



A small deviation from criticality, to oc A5 — Ap), intro- 
duces a term proportional to (to zeroth order in the 



small momentum k^), which identifies the remaining two 
gamma matrices. 



74 



75 



T a' 



(12) 



where 75 = 71727374 was used to determine the last 
gamma matrix. It is now indeed clear that 75 is odd 
under both time reversal (because it is proportional to 
spin (T^) and inversion (it is odd under layer exchange 
— > — T^), as argued on general grounds in the previous 
section. 

From that general analysis, we can view the critical 
Dirac state as the "parent" state of the topologically- 
stable nodal semimetal phases, which we will consider be- 
low. The simplest and most robust such phase is the Weyl 
semimetal, which we consider first. We saw in Sec. |IIBT] 
that the vector TR symmetry breaking perturbation b 
robustly splits the Dirac point into two Weyl points along 
the axis, parallel to b. Splitting the nodes along the z 
axis therefore is accomplished, from Table |lj by adding a 
term proportional to 6^ = 712 = , using Eq. ( [IT| ). This 
is precisely the spin-splitting term considered by two of 
us in Ref. |9l After a canonical transformation: 



cr — > T cr 



— cr r 



(13) 



the momentum-space Hamiltonian of the multilayer can 
be written in a block-diagonal form, with two indepen- 
dent 2x2 blocks: 

■H(k) = vpkya"" - vpk^a^ + m± {kz)a\ (14) 

where m±{kz) — b± A{kz), b is the coefhcient of the 
term (magnitude of the spin-slitting) and: 



A{kz) = ^jAi; + Ajy + 2AsADCos{kzd). (15) 

Taking 6 > 0, the to+ mass in always nonzero, corre- 
sponding to a pair of fully gapped bands. The to_ mass, 
on the other hand, changes sign at kz — n/dztko, where: 

fco = J arccos [l - (5^ - (As - Ac)2)/2AsAz5] . (16) 

The two points, where to_ vanishes, correspond to the 
two Weyl fermions, separated in momentum space. The 
Weyl semimetal phase exists as long as: 

bl, = {As - Anf <b'< bl, = (As + An)'. (17) 

As discussed in Ref. i9i Weyl semimetal is characterized 
a finite Hall conductivity, proportional to the separation 
between the Dirac nodes: 



(18) 



and chiral edge states, which exist only in a finite subset 
n/d — ko < kz < n/d + ko of the 2D BZ of any sample 
surface, not normal to the 2-axis. For more details on 
this we refer the reader to Ref. [HI 
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B. Effect of orbital field 

The realization of a Weyl semimetal, we have described 
above, requires doping the TI layers with magnetic im- 
purities. This is needed to produce the spin-splitting b, 
which breaks TR symmetry and splits the massive Dirac 
fermion into two massless Weyl fermions. In practice, 
it is easier to break TR by simply applying an external 
magnetic field, instead of doping the multilayer material 
with magnetic impurities. In this section we will explore 
this route in some detail. 

Let us assume that an external magnetic field of mag- 
nitude B is applied along the z-axis, i.e. the growth 
direction of the multilayer. The Hamiltonian is given by: 



Sct^-KA, (19) 



where 

A = A5t"<5.j + ^ (r+(5,- ,+1 + r-^,- , (20) 

is the tunneling operator in real space. We choose Lan- 
dau gauge for the vector potential A = xBy. Since the 
vector potential does not enter in the tunneling term A, 
it can still be partially diagonalized by Fourier transform. 
Then, after the canonical transformation of Eq. ( 13 ), and 



after diagonalizing the tunneling term, we obtain: 

H = vp{z X cr) • (^-iV + ^a) + m±(fcJcr^ (21) 

where m±{kz) = b± A(fcz), and b = gfiBB/2. This is 
identical to the problem of 2D Dirac fermions with masses 
m±{kz), which depend on a parameter fcz, in a perpen- 
dicular magnetic field. The solution for the spectrum is 
well-known and is given byi^it^il 



(22) 



where A = ± labels the electron and hole-like sets of 
Landau levels, ojb — vpjl is the analog of the cyclotron 
frequency for Dirac fermions, I = \fcfeB is the magnetic 
length (we will be using h = \ units throughout, restoring 
explicit Ti in some of the final results), and n — 1,2,... 
are nonnegative integers. As is well-known, the n = 0, 
i.e. the Lowest Landau Level (LLL), is special and needs 
to be considered separately. The energy of the LLL, cor- 
responding to the mass m+(fcz), which is always positive, 
is given by: 



eo+(fc2 



-m+(fcz), 



(23) 



i.e. the 0-1- level is always hole-like and lies below the zero 
energy line for any value of the momentum k^ . The sit- 
uation is different for the LLL of the Dirac fermion with 
the ■m-{kz) mass. m-{kz) = b — A{kz) changes sign from 
negative to positive as the momentum kz crosses the lo- 
cations of the Dirac nodes kz = n/d ± k^. This means 
that the 0— LLL is electron-like when jfe^ — 7r/o?| > fco. 



while it is hole-like, i.e. dips below the zero-energy line 
when \kz — T^/d\ < fco, or, in other words, when kz is in the 
interval between the Dirac nodes. This corresponds to a 
jump in the Hall conductivity of the corresponding ficti- 
tious system of 2D Dirac fermions, parametrized by fc^, 
from to /h. The total Hall conductivity of the multi- 
layer is obtained by integrating the 2D Hall conductivity 
between the Dirac nodes and is still given by the same ex- 
pression, as in the case of the magnetic-impurity-induced 
spin-splitting, without any orbital component of the field, 
i.e. (Jxy — e^fco/7r/i. The edge states also retain their 
character: these are chiral topologically-protected edge 
states, that exist in the interval ir/d—ko < kz < ir/d + ko 
in the 2D edge BZ. Thus some of the defining and most 
interesting properties of the Weyl semimetals can be ob- 
served by simply applying external magnetic field to a 
TI-NI multilayer structure, without any doping by mag- 
netic impurities. Also note that the magnetic field depen- 
dence of axy, which is given by Eq. (16), since axy ^ fcoj 



is quite different from what would be expected in a reg- 
ular metal. Indeed, ko is a highly nonlinear function of 
B, vanishing as ^/l)^^^JA^s^^^J^ near the transition 
from the Weyl semimetal to the insulator. 

Finally we remark that recent work has suggested 
a magneto-conductivity, i.e. diagonal conductivity for 
current and electric field parallel to an applied mag- 
netic field, for a Weyl semimetal. This was argued to 
be a manifestation of a "quantum anomaly" for Weyl 
fermions. While interesting, we note that a significant 
effect occurs only in the ultra-quantum limit in which 
ujbt ^ 1, where r is the scattering time. Moreover, a 
large conductivity, parallel to an applied magnetic field 
(relative to the orthogonal components) is in fact a rather 
generic consequence of the ultra-quantum limit, due to 
the quenching of kinetic energy in the transverse direc- 
tions and the suppression of backscattering in the effec- 
tively one-dimensional transport regime, resulting from 
high field. Thus association of this magnetoconductivity, 
parallel to the applied field, with Weyl physics, seems 
challenging experimentally. 



IV. LINE NODE SEMIMETALS 

In this section we will describe a realization, in the 
same physical system of a TI multilayer, of a line-node 
semimetal: a distinct topological semimetal phase, with 
zeros in the spectrum, forming continuous lines in mo- 
mentum space. 



A. Parallel-field-induced nodal line 

We consider a TI multilayer system in the presence of a 
magnetic field, parallel to the layers. This can be a real 
external magnetic field, or, as in the previous section, 
an exchange field, arising from ferromagnetic ordering 
of magnetic impurities, introduced into the TI material. 
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In the case of an externally applied field, we, for now, 
neglect the orbital effect of the field, but will discuss it 
in detail later. We anticipate the presence of nodal lines 
from the Dirac calculations in Sec. |IIB 2[ Indeed, from 
Table |l] and Eqs. ( 11 ), ( 12 ), an in-plane field corresponds 
to X and y components of the b' perturbation, which 
leads to a circular node in a plane containing the z axis. 

Let us consider this in more detail. Without loss 
of generality we assume the field is applied in the x- 
direction. The momentum-space Hamiltonian is given 
by: 

H(k)=WFr^(zxcr).k-H6a^-HA(/c^), (24) 

where 

= AsT- + i(AoT+e''^-^'^ + h.c). (25) 

The corresponding band dispersion, obtained by diago- 
nalizing Eq.(24), is given by: 



(26) 



where A{kz) is given by Eq. ( 15 ). The e_ branch exhibits 



a line node in the yz-plane, given by the solution of the 
equation: 



vlkl + A\k^)^b\ 



(27) 



As above, for concreteness we assume that As.d > 0. 
Then the node will be centered a,t ky — 0, k^ ^ n/d. The 
node exists as long as: 



b> |As- Aul. 



(28) 



B. Stability of the parallel-field-induced nodal line 

1. Nodal lines in general 

While we do not expect complete stability of the nodal 
line, it could be stabilized if extra symmetries are im- 
posed upon the Hamiltonian. We need to distinguish 
two types of stability. First, we can ask whether the line 
contact of conduction and valence bands is stable. Sec- 
ond, we can ask whether, if this is stable, the line contact 
is degenerate and coincides with the Fermi energy. The 
answer will be that the former is possible with some dis- 
crete symmetries, while the latter cannot be guaranteed 
by any set of discrete symmetries, although it will be ap- 
proximately degenerate, and perhaps to a high degree of 
precision, under most reasonable circumstances. 

First we discuss the question of band touching from a 
general point of view. Since the nodal line occurs in a 
system with non-degenerate bands (away from the node 
itself), it is sufficient to consider a two-band Hamiltonian. 
This in general takes the form: 



where the Pauli matrices cr act in the two-band space. 
By simple counting of the degrees of freedom and con- 
straints, nodal lines may occur when, for instance, one 
of the ha, for a ~ 1,2,3, vanishes for all k. Then as a 
function of k, two parameters must be tuned to obtain 
band degeneracy, i.e. to make the other two ha vanish, 
which results in line nodes in momentum space. Note 
however, that even in this case, any momentum depen- 
dence of /io(k), which in general is not constant along the 
contact line, means that the line touching does not have 
a constant energy, and therefore cannot coincide with the 
Fermi level. In general, there are an infinite number of 
functions /io(k), consistent with any discrete symmetries, 
so that this requires an "infinite" degree of fine-tuning. 
This means that a line-node semimetal, i.e. a semimetal 
with a line-like Fermi "surface" in 3D is non-generic. 

Nevertheless, line contacts, even with non-constant en- 
ergy, have a robustness associated with them, which can 
be traced to the existence of a topological invariant, char- 
acterizing the line contact. For concreteness, and with- 
out loss of generality, let us take hi ~ 0. Then we may 
form a complex order parameter h — h2 + ih^, whose 
phase is well-defined everywhere except at a node. Away 
from the node, we define h = |/i|e*^, and then, since h is 
single-valued, we have: 



(30) 



where n is an integer for any closed curve C in momentum 
space, on which the bands are non-degenerate. Since this 
winding number is quantized, it cannot change as this 
curve is smoothly deformed. If the curve does not contain 
any singularity inside it (i.e. points where h vanishes - 
nodes), then it can be shrunk to a point and the winding 
number n must vanish. Generically, however, a curve 
that encircles a nodal line, has n — ±1, depending upon 
the sense of circulation. The nodal line can therefore be 
viewed as a vortex line in momentum space. So long as 
vanishing h requires a band degeneracy, this is the case. 
However, if we allow the third component, hi ^ 0, then 
h = h2 = h^ = does not require a band degeneracy, 
and a curve with non-zero n need not enclose a node. 

This can be understood in a yet more general context. 
Specifically, the line integral in Eq. (30) can be viewed 
more generally as a Berry phase. Whenever the bands 
are non-degenerate, we can define a U{1) Berry gauge 
field (Berry connection) from the periodic part of the 
Bloch wavefunctions UkQ(r), where a is the "spin index", 
associated with the Pauli matrices in Eq.(29): 



dk, 



r^u.c. 



d 



-wL(i")wka(r) 



(31) 



'H-2b{k.) — /io(k) -I- /ii(k)CT^ + /i2(k)(T^ + /i3(k)(T^, (29) where the integral is taken over the unit cell of the crystal. 
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The Berry curvature is the flux of this gauge field: 

-B^(k) = e^,,xd,Ax (32) 



Then we obtain two independent blocks of the Hamilto- 



d 



»7^Wka(r) 



Stable nodal lines occur when the Berry curvature is 
generically (i.e. for non- degenerate points) vanishing. 
This is because one may write the line integral as: 



i 



(33) 



by Stokes' theorem, where iS is a surface in reciprocal 
space, whose boundary is C. For any surface, for which 
there is no band touching, = would imply a van- 
ishing "vorticity" . Conversely, non- vanishing vorticity 
within C implies non-vanishing Berry curvature on S. If 
the Berry curvature is generically zero, then this in turn 
requires a singularity on S, i.e. that S is crossed by a 
nodal line. However, if there is no such requirement of 
vanishing Berry curvature, there need be no singularity, 
and the curvature may be spread out over the region of 
integration. 

The vanishing Berry curvature condition holds in the 
above example because, when hi — 0, the Hamilto- 
nian obeys a^H* (k)a^ = H(k). When this condition is 
obeyed, the Bloch functions satisfy = cr^^Ukp, which 
implies a vanishing Bf^ from Eq. ( 32 ) . One can see that 



such a vanishing-Berry-curvature condition generally re- 
quires some discrete symmetry, involving conjugation of 
the Hamiltonian at a single momentum point. Without 
both inversion and time reversal symmetry present, this 
is, in general, artificial. Nevertheless, it may be imposed 
in toy models, or may be approximately the case for some 
physical situations, such as discussed here. 



2. Superlattice case 

We now return to the specific case of the nodal line, 
induced in the TI-NI superlattice by an in-plane field. It 
is instructive to reduce the Hamiltonian to a two-band 
form, containing just the bands involved in the line node. 
To do so, we first rotate the spin quantization axis by 7r/2 
around the y-axis, taking % ^ %: 

H{k) ^{b + vpT'ky) a' + A(fc,) - VFT'<jyk^, (34) 



and then make the canonical transformation of Eq. ( 13 1, 
under which "H — >■ "H', with: 



n'ik)= b + VFT'ky + A{k,) 



(35) 



The term in the square brackets is now a constant of 
motion and can be replaced by its eigenvalues: 



'H'-i-(k) — m±a^ — ypk^cr^. 



(37) 



The low-energy block, containing the node, corresponds 



to 'H'_, and indeed has the form, described in Sec. IV B 1 
It has a symmetry W (for "wishful thinking" ) : 



W : a' [H'±(k)]*(7" =Hi(k). 



(38) 



But is this symmetry physical? 

Complex conjugation occurs physically only through 
time reversal, which we denote by T. T acts on the 
original Hamiltonian, Eq. (|24]), as 'H(k) -> cr^'H* (-k)cr2'. 
It is not a symmetry due to the applied Zeeman field h. 
However, if combined with a tt rotation about the z axis, 
7?.|, the invariance is restored. So a physical symmetry 
is: 

ToUl: H{k,,ky,k,) ^ (T^H*{k^,ky,~k,)a\ (39) 

Carrying through the transformations from 'H(k) to 
'Hj-(k), invariance under To TZ^ requires: 

ToUl : a' [n'±{k,,ky,k,)Y =n'±{k,,ky,-k,). 

(40) 

This is close to, but not precisely the required condition 
to protect the node, because it involves a sign change 
of kz- We can attempt to reverse this sign change by 
imposing an additional z — ^ — z reflection symmetry on 
the problem, which might naturally be associated with 
reflection with respect to a constant z plane at the cen- 
ter of the TI or NI layer. However, because spin is a 
pseudovector, this will flip the in-plane components of 
the spin tr^ — ?► — c^, cr^ — ^ —a^. This does not leave the 
in-plane field ba^ invariant, so cannot be a symmetry of 



Eq.(24) 



Because the growth direction and applied field fully 
break any possible three or fourfold rotation axes, only 
some discrete Z2 type symmetries remain as candidates. 
One possible remaining symmetry consistent with the 
applied field is inversion, I, through a center mid-way 
through a TI or NI layer. Since spin is a pseudovector, 
this leaves cr invariant. Upper and lower layers of each 
TI layer are interchanged, and k — >■ — k, so this condition 
gives: 



X : nik) ^ T^-H(-k)T^, 



(41) 



which is indeed an invariance of Eq.(|24[). After the 
changes of basis, I implies for Eq.(35l that: 



(42) 



Finally, after projection into the 2x2 blocks, the final 
condition, imposed by inversion symmetry I, is: 



m±(k) = 6±./4/c2 + A2(fc, 



(36) 



I: CT"H'±(k)(7" =H'±(-k). 



(43) 
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Another possible symmetry is a two-fold (tt) rotation 
about the axis of the field, TZ^. This acts as: 

: n{k.,, ky, k,) ^ a^T-H(fc„ -ky, -k,)T-a\ (44) 



number analogously to Eq. ( 30 ) , but with the contour C 



which can be verified to be a symmetry of Eq. ( 24 ) . Map- 
ping this to the rotated frame, we obtain: 

: n'ik,, ky, fcj ^ T-n\k,, ^ky, ^k,)T\ (45) 

Projecting down to the 2x2 blocks, this gives: 



'^ir • 'H±{kx , ky , kz) — Hj-i^kx, ky, kz 



(46) 



Even if invariance under all three symmetries (To72.|, 
I, and TZ%) is imposed, this is not equivalent to Eq. (38 1, 
and the nodal line is, at first sight, not guaranteed to 
be stable. In particular, a non-zero /ii(k) (proportional 
to CT^) is allowed to enter in Eq. ( [37| , which could 
destabilize the nodal line. However, one may show that 
/ii(k) must, according to these three symmetries, be an 
odd function, separately, of k^, ky and kz- In particular, 
this implies that the hi{kx = 0, ky, kz) = term vanishes 
on the y~z plane, where the nodal line exists. Therefore, 
the hand contact along the nodal line is indeed protected 
when all three symmetries are present. 

We may also consider, however, the constant part of 
the two-band Hamiltonian, /io(k). This is required, by 
the same symmetries, to be an even function, separately 
of kx,ky, and kz, and hence does not vanish nor need be 
a constant at kx = 0. Physically, this term would arise, 
for example, from the always-present particle-hole asym- 
metry of the TI surface states. Generically this splits the 
zero energy Fermi line into a set of small electron and 
hole Fermi surfaces, converting the line node state into 
a conventional low carrier density semimetal. However, 
the topological surface state, associated with the nodal 
line, survives the addition of ho(k.) term, although ac- 
quires a dispersion (while it is strictly dispersionless in 
the absence of /io(k)). 



C. Surface states 

While the line node at the Fermi energy is not generic, 
it may well be a good approximation, and regardless, the 
line band contact itself is more robust, as we have seen 
above . It is interesting to understand the consequences 
of this bulk band topology for the boundary. Indeed, 
Heikkila and Volovik have shown in another context that 
un usual surface states are related to a nodal degener- 
g^pylllSOl uncover the nontrivial surface effects of the 
line node, it is useful to view the Hamiltonian 'H(k) as de- 
scribing a set of ID systems, parametrized by momentum 
components ky,kz'. Hk ,kA^x)- Such a one dimensional 
two-band Hamiltonian supports a topological classifica- 
tion, if only two of three Pauli matrices are present in 
it, which is the same as the condition to generically sup- 
port a line node. In this case, we can define a winding 



taken along the periodic direction kx in reciprocal space, 
i.e. across the entire BZ. Thus the same condition, which 
generically gives stable nodal lines, also allows such a ID 
topological classification. When the ID winding number 
n is non-zero, then a bound state is expected at an in- 
terface between the system and another system with a 
different value of n, e.g. n = 0, corresponding to the 
vacuum. Indeed, recall the formula for the canonically 



transformed Hamiltonian in Eq. (37 1, but regarded as a 



ID Hamiltonian, parametrized by ky,kz'. 

n'_.k,,kAkx) = m.iky,kz)a'- - vpaykx. (47) 

The mass TO_(k) changes sign from negative to positive 
when: 



b,{ky,kz) = ^Vj,kl+A^kz). 



(48) 



Although this Hamiltonian is written only for small kx, 
and hence does not describe the full ID topology, it does 
describe transitions between different topologies, which 
occur when the mass m_ changes sign. The ID TI-NI 
transition occurs when the above condition is satisfied, 
so that when m_ > 0, one has a non-trivial ID insulator 
and surface bound states, while when m_ < the ID 
insulator is trivial and no bound states are guaranteed 
at the surface for such ky,kz. 

To see this explicitly, let us assume that the sample 
occupies the a; < half-plane with a surface at a; = 0. To 
find the edge states of the ID TI inside the nodal line, we 
replace kx —id/dx and look for solutions of the 2x2 
Dirac equation: 



•H' .j^ fe * = 0, 



in the following form: 



(49) 



(50) 



where is a two-component spinor. Substituting this 
ansatz into the Dirac equation, we obtain: 



TO_ {ky, kz,x)a^ — vp 



dF 
dx 



= 0. 



(51) 



Assuming b{x — ^ oo) = 0, the solution is given by: 

^k,,kM = ,^i,ro'^-'^-i^y,^^:-\. ^ (52) 

This is normalizable and localized at the surface for all 
bc{ky,kz) < b. The set of these zero-energy edge states 
forms a flat band in the surface BZ, which is dispersion- 
less for all ky, kz inside the area, enclosed by the projec- 
tion of the nodal line on the surface plane. If we now 



add the Hq (k) perturbation to the Hamiltonian Eq. ( 47 1 , 
it is easy to show, using the standard quantum mechan- 
ical perturbation theory, that to leading order in Hq, 
the surface state acquires a dispersion, proportional to 
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FIG. 1. The fcz = 57r/6d section of the eigenstate spectrum for 
a sample of finite size in the a;-direction with Ad/As — 0.8, 
and b/As ~ 1. ky is in units of As/vp. The intensity of gray 
is a function of the degree of surface localization of a given 
eigenstate, measured by an inverse participation ratio of its 
wavefunction. The surface state dispersion is black, while 
bulk states are lighter gray. 




FIG. 2. The kz = 57r/6d section of the eigenstate spectrum for 
a sample of finite size in the a;-direction with Ad/ As ~ 0.8, 
and b/ As = 1 in the presence of a particle-hole asymmetry in 
the TI surface state spectrum of the form {k^ + k^) /2m' , with 
1/m* — Q.?iv% / As- ky is in units of As/vp- The surface state 
has acquired a dispersion due to the particle-hole asymmetry. 



hQ{0,ky,kz). The full surface fiat band dispersion, cal- 
culated numerically (in the absence of the ho term), is 
shown in Fig. [l] 

In the presence of /io(k) terms, that will give the sur- 
face state a dispersion, the remaining robust topological 
property of the surface state will be its termination at the 
projection of the nodal line to the surface BZ, as shown 
in Fig.[2| 

As discussed above, ho(k) kills the bulk line node itself, 
transforming the line-node semimetal to a more "conven- 
tional" semimetal with a Fermi surface, containing elec- 
tron and hole pockets of equal volume (at charge neu- 
trality), as shown in Fig. [s] It is, however, distinguished 
from a truly conventional semimetal by the presence of 
the topological surface states, described above. 




FIG. 3. (Color online). The "Fermi line" transforms into a 
finite-volume Fermi surface with equal-volume electron (top 
and bottom) and hole (left and right) pockets due to the 
particle-hole asymmetry of the TI surface states. Note that 
the line contact of the conduction and valence bands survives, 
though it is not a constant energy curve. This line contact 
threads through the middle of the chain of small Fermi sur- 
faces, and as a consequence the topological surface state re- 
mains (see Fig. [2 1. The parameters here are taken to be the 
same as in Fig. 



D. Effect of the orbital part of the field 

We have so far ignored completely the effect of the or- 
bital part of the parallel field, which creates the line node 
in our system. This is justified when the parallel field is 
an exchange spin-splitting field, coming from the interac- 
tion with ferromagnetically-ordered magnetic impurities. 
If the field is an externally-applied magnetic field, how- 
ever, a situation that is perhaps more easily realizable 
experimentally, the orbital effect of the field needs to be 
considered, as in the case of the Weyl semimetal, dis- 
cussed above. 

The Hamiltonian in the presence of a magnetic field of 
magnitude B, directed along the x-axis, is given by: 



H ~ vpT^iz X cr) 



Ba^' + A, (53) 



where A is the tunneling operator in real space, given 
by Eq. ( [20| ). We choose Landau gauge for the vector 
potential A = —zBy, in which case A does not enter 
in the tunneling term. We would like to point out here 
that the magnitude of the ij-factor is large in typical TI 
materials, e.g. t; « 50 in Bi2Se3, so the Zeeman term will 
produce a significant spin-splitting at reasonable values 
of the magnetic field. This remark is also relevant for the 
Weyl semimetal in magnetic field case, discussed above. 



10 



The general case of TI multilayer in a parallel field can 
only be studied numerically, as this is a Hofstadter-type 
problem. However, two limits can be studied analytically. 



1. Limit of almost decoupled TI layers: Ad ^ As 

In this limit the problem reduces to the one of a single 
TI layer in parallel magnetic field, already considered by 
two of us in Ref. Neglecting the contribution of A^i, 
A(fcz) = As becomes independent of momentum. The 
nodal line in this case has the form of two straight lines, 
parallel to the z-axis, crossing the y-axis at ky — ±fco, 
where: 



(54) 



cb — vpUB '^^ the "magnetic energy", with the "magnetic 
wavevector" kb given by: 



(55) 



where dri is the thickness of a TI layer and £ — \J cjeB 
is the magnetic length. The first term in Eq. ( 55 1 comes 



from the orbital part of the field, while the second term 
comes from the Zeeman spin-splitting part. The topolog- 
ical surface state in this case consists of a set of ID edge 
states of each TI layer, which are dispersionless in the 
y-direction when —ko<ky<kQ. The surface state also 
does not disperse in the z-direction, but for the trivial 
reason of the absence of tunneling between individual TI 
layers. 



2. Weak field limit 

The weak field limit applies when: 

\As-Ad\^As + Ad, (56) 

which implies that we can assume: 

6<As + Ad, £>d, (57) 

where we have used short-hand notation for the Zeeman 
spin-splitting term: 



(58) 



In this case, first setting the orbital part of the field to 
zero, we can expand 'H(k;) in Taylor series around k^ = 
Tr/d, which is the location of the nodal line center in the 
absence of the orbital component of the field. Expanding 
to leading order in k^ — Tr/d, and shifting the zero of the 
momentum to k^ ~ 7r/d, we obtain: 



A(fc,) « (As - Ad) + Andryk, 



(59) 



Rotating by 7r/2 around the y-axis and performing the 
canonical transformation of Eq. ( 13 1 , the problem reduces 



to finding the Landau level spectrum of the following 2x2 
Hamiltonian: 

H = vpT'TTy + Acdr%^ -I- (As - A^) t\ (60) 
where tt is the kinetic momentum: 



._9 _ ^ 



._d_ 
dz 



(61) 



The Landau level spectrum is easily found in the stan- 
dard way, by introducing ladder operators as: 



VF 

2vpP 



{a} + a) 



VF 



(at 



(62) 



where vp — d^/AsAjj ^ dAjj has the meaning of the 
Fermi velocity, associated with the z-direction in momen- 
tum space. The resulting Landau level spectrum, taking 
As ~ Ad, is given by: 



±A/2w|n, 



(63) 



where ujb — \JvpVp j i and n = 0, 1, 2, . . . The full Hamil- 
tonian Eq.(53l can now be written as: 



b± X 2iu 



— vpu'^k^. 



(64) 



By exactly the same reasoning as in section |IV A[ we 
can conclude that topological zero-energy surface states 
appear when: 



b > X 2u 



(65) 



The surface states in the case of an externally-applied 
parallel magnetic field will thus consist of Landau levels, 
which become localized at the surface of the sample, nor- 
mal to the applied field. In other words, the topological 
surface flat bands at zero field transform into surface- 
bound Landau levels in an applied magnetic field. 



V. CONDUCTIVITY OF THE NODAL STATES 

In this section, we discuss the DC and optical con- 
ductivity of the nodal states described in the previous 
sections. 



A. Weyl semimetal 

Here we focus on the diagonal transport characteris- 
tics of the Weyl semimetal, namely its optical conductiv- 
ity. The simplest possible calculation of the conductivity. 
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neglecting interactions, assuming charge neutrality, and 
taking into account only random point impurities, was 
quoted in Ref. IHl but the results were not derived in de- 
tail. Here we give a somewhat more general discussion, 
and in particular show that, in fact. Coulomb interac- 
tions drastically change the behavior of the conductivity 
at low temperature. In particular, the result of Ref. P 
survives only at high temperature, and if the basic inter- 
action scale, defined by the effective fine structure con- 
stant, a = /edVp, where is the dielectric constant, is 
small, a ^ 1. In general, we will show that the frequency, 
doping and temperature dependences of the conductiv- 
ity of the Weyl semimetal are very unusual, and can be 
used for experimental characterization of this phase of 
matter. Moreover, with interactions and charged donors 
taken into account, we argue that the conductivity obeys, 
up to logarithmic corrections, a scaling form in its depen- 
dence upon T, uj, and donor impurity density n^. 

The failure of the non-interacting point impurity re- 
sult is in dramatic contrast to conventional metallic sys- 
tems, in which elastic scattering from defects dominates 
over inelastic electron- electron processes at low tempera- 
ture, which are frozen out due to phase space restrictions. 
Even in the apparently close analog of 2D graphene, dis- 
order dominates the low energy transport rather than in- 
teractions, in striking contrast to the 3D Weyl semimetal. 
This difference can be explained without detailed calcu- 
lations from a simple renormalization group (RG) argu- 
ment. Consider the action of the Dirac/Weyl fermion 
model in d dimensions with both point disorder and in- 
teractions with a 1/r Coulomb potential: 



drd^x 



■0(9r + i^t,d^)-ip + V^{x)'ipAiiJj 



+ / dTd^xd''x'{4,^P),^r- (66) 

J 2\x-x'\ 

where here Vi{x) are random potentials, coupling to the 
fermion fields via some matrices Ai (not specified), and 
e is the electron charge. We take the quenched ran- 
dom potentials to have zero mean and Gaussian vari- 
ance Vi{x)Vj{x') = AijS^'^\x — x'), reflecting short-range 
correlations. To keep the free Dirac/Weyl action scale 
invariant, we must under an RG transformation rescale 
length and time as a; — )■ bx, t ^ bt and ip b~'^^'^il}. Un- 
der this rescaling, we see that the Coulomb interaction 
term, proportional to e^, is marginal in any dimension. 
However, the random potential Vi{x) — ^ bVi{bx)^ which 
implies that the disorder strength A^- — ^ b'^^'^/^ij. Thus 
in the d = 2 case of graphene, disorder and interactions 
are both marginal by power counting. In fact, more care- 
ful analysis shows that interactions are marginally irrele- 
vant and disorder is marginally relevant at the free Dirac 
fixed point. As a consequence of the marginally relevant 
disorder, a density of states is generated and the sys- 
tem is described at low energy by a diffusive fixed point. 
By contrast, for the d = 3 Dirac/Weyl fermion, interac- 
tions remain marginal (actually marginally irrelevant, as 
a more detailed analysis showP^), but disorder becomes 



strongly irrelevant A — > A/6. Thus in fact the ballistic 
fixed point is stable for weak disorder in three dimen- 
sions. Moreover, since disorder is much more irrelevant 
than interactions at this fixed point, elastic scattering 
is suppressed relative to inelastic scattering, and this ex- 
plains the dramatic difference from the 2D graphene case. 
In fact, the marginality of interactions means that many 
physical properties almost scale like those, expected of a 
fully interacting scale-invariant critical theory, with only 
logarithmic corrections. This simple "quantum critical" 
scaling is an attractive feature of the Weyl semimetal. 

In the remainder of this subsection, we will go beyond 
these scaling considerations and verify their conclusions 
in some simple calculations. We will also extend the 
discussion to the physically relevant situation, in which 
donor impurities are present, which extends in a simple 
way the quantum critical scaling due to interactions to 
include finite residual resistivity at T = 0. 



1. Short-range impurities for non-interacting electrons 

We first recapitulate the calculation of Ref. [9l since 
all details were omitted in it. We assume a model with 
short-range impurity scattering potential of the form: 



(67) 



where label the impurity positions, and also neglect 
electron-electron interactions. Both of these assumptions 
are generally quite unrealistic, both for an undoped Weyl 
semimetal, in which Coulomb interactions are essentially 
unscreened and for a doped semimetal, where scattering 
from charged donors with long-range potential can be 
expected to dominate. However, this model will still give 
us useful results, which can be expected to be applicable 
at the neutrality point at high enough temperature, such 
that the impurity scattering rate exceeds the scattering 
rate due to electron-electron interactions. 

We will assume that the impurity potential is diag- 
onal in both the spin and the pseudospin indices and 
will consider a single Weyl fermion in the 3D BZ with a 
Hamiltonian: 



?^(k) = vpcT ■ k. 



(68) 



Generalization to any number of distinct Weyl fermions is 
trivial, as they contribute additively to transport (we will 
assume that the impurity potential does not mix Weyl 
fermions at different points in the BZ). 

In the first Born approximation, the impurity scatter- 
ing rate is given by: 



1 



r(6) 



where 



= -7 J 



/^3l ^ ^ 
^^Gf(e,k) = 27r7g(6), 



1 



e — Xvpk + irj 



(69) 



(70) 
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is the retarded Green's function of the Weyl fermion, 
A = ± labels the helicity of the positive and negative 
energy Dirac cones, 7 — u^rii, where is the impurity 
concentration, and the density of states g{t) is given by: 



5(e) 



(71) 



Thus 1/T(e) ^ <C e, which means that the conduc- 
tivity can be calculated semiclassically, using Boltzmann 
equation. Solving linearized Boltzmann equation with 
the energy-dependent momentum relaxation rate ( 69 1 in 
the standard way, we obtain: 



Re axx{uj) 



de g{e 



dnpie) l/r(e) 



(72) 

where np is the Fermi distribution function at temper- 
ature T. Introducing dimensionless integration variable 
X — e/2T (using fc^ = 1 units) and restoring explicit Ti, 
we obtain: 



Re axx{(^) 



""""P I dx 



67/1 



X* sech^(a:) 



x4-|-(/i3i;3cj/327r27T2)2' 



This gives a DC conductivity: 



ctdc 



C Up 

3jh ' 



(73) 



(74) 



and a Drude-like peak in the optical conductivity, but 
with a temperature-dependent width, scaling as . This 
is a very unusual property of the optical conductivity 
in a metal and can be used to characterize the Weyl 
semimetal phase experimentally. 

The Drude peak also has a highly unusual shape, with 
a divergent first derivative. This can be obtained explic- 
itly from Eq. (73), evaluating the integral in (731 in the 
limit w — > 0: 



Re a^x{LLi) 




(75) 



2. Donor impurities 

Now let us consider a doped Weyl semimetal and 
adopt a more realistic model with Coulomb, rather than 
short-range impurities, which will represent the charged 
donors. As is well-known, Boltzmann approach can still 
be used in this case, with transport time replacing the 
momentum relaxation time of Eq.(69): 



-—- =7rn,g{e) / dO sini0)\V iq)\^ 

Ttr[e) Jo 



X [1 - cos(6')] 



1 + cos(0) 



(76) 



where 



(77) 



is the screened Coulomb potential with the Thomas- 
Fermi wavevector q^p — 4:7Te'^g{e), 9 is the scattering 
angle, q — 2(e/wi?) sin(0/2), and ni is the impurity con- 
centration. The factor 1 — cos(6') is the standard factor, 
suppressing the forward-scattering contribution to the 
transport collision rate, while the (1 + cos((?))/2 factor 
arises from the matrix elements of the impurity poten- 
tial with respect to the eigenstates of the Weyl Hamil- 
tonian Eq.(68|. Eq.(76) is in fact very similar to the 



corresponding expression for the transport collision rate 

in graphene.'^SJ 

Introducing an effective fine structure constant a = 
je^vp^ which expresses the ratio of the typical Coulomb 

interaction energy scale e^kp/cd to the typical kinetic 

energy scale vpkp in the Weyl semimetal, Eq.(76l may 

be written as: 



1 



Ttr(e) 



■na^riiv'p 
4e2 



dO 

[sin2(6'/2) +a/27rl 



(78) 



Integrating over the scattering angle, we then obtain the 
following expression for the transport collision rate: 



1 



47r^r 



'J '^p 



3e2 



where 



/(«) 



3a2 



[I 



i!/7r)atanh 



1 + a/TT 



1 



(79) 



(80) 



The function f{a) approaches unity for a 3> 1, i.e. in 
the limit of strong interactions and vanishes as: 



30-2 

/(a)«^ln(l/a). 



(81) 



in the weak interaction a ^ 1 limit. If we assume that 
the charged impurities are donors, i.e. that the impurity 
concentration is proportional to the electron concentra- 
tion n ^ {ep/vp)^, then, as obvious from Eq.(79l, 



/(«) 



1 



(82) 



The standard Boltzmann-cquation expression for the 
zero-temperature DC conductivity in terms of the trans- 
port collision time is given by (h is restored): 



9 2 

r^DC ~ —r^9{>^F)Ttr{e-F)- 



(83) 



Using Eq.(79) and Ui {ep/vpy, we finally obtain the 



following result for the DC conductivity of the Weyl 
semimetal with Coulomb impurities: 



2 1/3 



(84) 
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Scattering from charged impurities thus leads to the DC 
conductivity vanishing as a function of the dopant con- 

1 /3 

centration as n - ' . Note that this result can be rewritten 
in the physically transparent form: 



(85) 



1 /3 

where the Fermi momentum kp ~ and t — vpT is 
the mean free path. This agrees with a simple scaling of 
conductivity linearly with energy. Moreover, the mean 
free path obeys: 



kpi^l/f{a) 



(86) 



which implies that electrons are weakly scattered and 
justifies the semiclassical approximation when /(a) <C 
1, i.e. when a <C 1, consistent with the perturbative 
treatment of scattering. 



3. Coulomb scattering at neutrality 

This behavior, however, can not be expected to hold 
all the way to the neutrality point, and a crossover to 
a doping-independent value of the conductivity should 
occur once ep < T, as happens e.g. in graphene.^^ In 
the regime ep < T, the conductivity can be expected 
to be determined by scattering due to the (almost) un- 
screened Coulomb electron-electron interactions (except 
at higher temperatures, where the short-range scattering 
from neutral defects, discussed above, will dominate). 
For undoped Weyl semimetal at low temperatures one 
expects the electron self-energy due to interactions to be 
proportional to the quasiparticie energy, up to possible 
multiplicative logarithmic corrections?^ 



1/t - Imi;(e) - a^Max{e,T}. 



(87) 



This follows simply from the absence of any energy scales 
in the undoped Weyl semimetal, other than the e itself, 
but can also be obtained from an explicit calculation.!^ 
Interpreting ImS(e) as the scattering rate 1/T{e) and 
plugging it into the Boltzmann-equation expression for 
the DC conductivity, we obtain: 



0) 



9 2 



de 



dnp{e) 



Vp 



i.e. a power-law insulating behavior with the DC con- 
ductivity vanishing linearly with temperature. 

Note that this form matches nicely the scaling, ob- 
tained above for the case of donor impurities. In partic- 
ular, both forms and the expected frequency dependence 
can be encompassed by the general quantum critical scal- 
ing form: 



a{ui,n^,T) 



e^kp 



S[vpkp/T,uj/T] , 



(89) 



where S[X^ Y] is an 0(1) scaling function. We have taken 
the forms appropriate for small a, where the perturbative 
calculations are valid, and neglected logarithms (which 
are interesting, but beyond the scope of this work). A 
more detailed study of the conductivity in the absence of 
donor impurities can be found in a recent preprint 



B. DC conductivity of the line-node semimetal 

Finally, let us discuss transport properties of the line 
node semimetal. As will be shown below, these are some- 
what similar to graphenej^l^^ except for the fact that 
they occur in a 3D material. We will also only consider 
the DC conductivity, as the calculation of the frequency- 
dependent conductivity is somewhat complicated and in 
general can only be done numerically (except in the high- 
frequency limit UJT ^ 1). However, as in graphene, we ex- 
pect the optical conductivity of the line-node semimetal 
to be only weakly frequency-dependent. 

For simplicity, we will adopt the model of point-like 
randomly-distributed impurities, with the potential given 
by Eq.(67). We assume that only the low-energy states, 
described by the lower (-) block of the Hamiltonian 



Eq.(37) contribute significantly to transport. The eigen- 
states of this 2x2 Hamiltonian are given by: 



±,k) = -^(Vl±™(k)/e(k), 

T isign(A:^)v/lT'7i(k)/e(k)) , 



(90) 



where m(k) = ?7i_(k) = 6 — ^ v^pk"^ + /S?[kz) and 



e(k) = Vpk"^ + m2(k). The corresponding eigenvalues 
are e±(k) = ±e(k). Let us first find the low-energy den- 
sity of states of the nodal line. In general it can only be 
calculated numerically. To obtain an analytical expres- 
sion we will assume that the size of the nodal line along 
the z-axis is small compared to n/d. Then we can expand 
A{kz) to leading order in k^ near k^ — n/d. We obtain: 



A(fc,) 



(91) 



where A = IA5 — Api\ and we have redefined k^ — > 
kz +Tr/d. The density of states can be most conveniently 
found by differentiating the function: 



7V(e) 



d^k 



3e(e-e(k)), 



(92) 



where Q{x) is the Heaviside theta-function. This is pro- 
portional to the volume of a torus in momentum space, 
whose surface is described by the equation: 



2 1,2 



,2 L2 



(93) 



The cross-section of this torus is not circular and its vol- 
ume in general can not be calculated analytically. An 
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analytical expression can, however, be obtained in the 
limit 6 ^ A. In this limit Eq.(93), after appropriate 



rescaling of the coordinates, describes a canonical torus 
of major radius b and minor radius e. Then we obtain: 



N{e) 



4:7TVpVF 



The density of states is thus given by: 



die) 



eb 



de 2TrvpVF 



(94) 



(95) 



The density of states of a 3D nodal line is thus the same 
(i.e. scales linearly with energy at low energies) as the 
density of states of point nodes in 2D, as expected. This 
means that many of the transport properties of the 3D 
line-node semimetal will be similar to those of graphene. 
In particular, since the first Born approximation scatter- 
ing rate l/r^e) ~ g{e) e, i.e. is of the same order as 
the quasiparticle energy, the first Born approximation is 
in fact inapplicable and the self-consistent Born approx- 
imation (SCBA) must be used instead. A general SCBA 
expression for the disorder self-energy is given gy: 



SA(k,e) 



where 



k',A' 



(VAA'(k-k')VA'A(k'-k))Gf,(k',e), (96) 



Gf(k,e) 



1 



(97) 



e-eA(k)-SA(k,e)' 
is the retarded disorder-averaged Green's function and 



VAA'(k-k') = (A,k|A',k')y(k-k'), 



(98) 



is the matrix element of the impurity potential. The 



angular brackets in Eq.(96) denote impurity averaging. 
Near the nodal line we can approximately set m(k) « 0, 
the matrix element in Eq. ( 96 1 becomes independent of 
A, A', and we obtain: 



E(6) 



2V 



EGA(k,e), 



(99) 



k.A 



i.e. the self-energy is independent of k and A. Since we 
are interested in S(e) at low energies, we can set e — > 
in Eq.([99|). Then we obtain: 



£2' 



(100) 



where S = E(0). Since g(e)is an even function of the 
energy, it follows from Eq. ( 68 ) that E is imaginary. Sub- 



stituting Eq.(|95[) in ([lOOj) and solving the resulting equa- 
tion, we obtain: 



1 

27 



-2ttv pVp 1 76 



(101) 



where Cc is the upper cutoff energy, which is of the order 
of the total bandwidth. Thus we find that the impurity 
scattering rate is finite in the zero-energy limit, unlike the 
naive first Born approximation result. It can, however, 
be very small in a clean multilayer. 

Now we can evaluate the conductivity. The standard 
Kubo formula expression reads: 



dt 



dnp{e) 



d^k 



^ J-00 V de J J (2^)3 
X (Ak|w„|A'k)(A'k|w^|Ak)ImGj?(k,e)ImGf,(k,e), 

(102) 

where repeated A, A' indices are summed over. Vertex 



corrections to Eq. (102) vanish identically. This can be 



checked by an explicit calculation, but is most easily seen 
Using the 2x2 momentum-space Hamiltonian: 



the components of the velocity operator v^, 
are given by: 



(103) 
dn/dk^ 



-vpa" 



Uprvy 



b — m(k) 



where at low energies we can again use m(k) 
elements of the spin operators are given by: 



6 — m(k) 

(104) 
0. Matrix 



-,k|a^ 
-,k|a^ 
-,k|a^ 



-,k) 
-,k) 
',k) 
-,k) 
-,k) 
^,k) 



e(k)' 

mCk^kx 

^(kT 
TO(k) 

m(k) 

VF\kx\ 

e(k) ■ 



(105) 



It follows from Eqs.(104| and (105) that at low energies, 
only intraband terms in (102) contribute to a^x, while 
only interband terms contribute to ayy and azz- After a 
straightforward calculation we obtain: 



nvph 



e^vpb 

TTVph ' 



(106) 



where we have restored explicit h. We note here that the 
vertex corrections to Eq. (102) vanish identically. This 



can be checked by an explicit calculation, but is most eas- 
ily seen from the following symmetry of the Hamiltonian 



Eq. ( |103[ ): -H(k) = n*i-k). It has been shown in Ref.[36] 
that such a symmetry of the Hamiltonian (combined in 
our case with the reality of all the matrix elements of 
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the velocity operator) always leads to cancellation of the 
vertex corrections to conductivity. 

The conductivity of the nodal-line semimetal is thus 
independent of disorder. This is similar to the well- 
known universality property of the DC conductivity of 
graphene.'^^^^ Unlike in graphene, however, the conduc- 
tivity of the line-node semimetal does depend on nonuni- 
versal properties of the nodal line, like its perimeter, 
which is proportional to b, and the Fermi velocity. Note 
that a XX = cTj/y only in the limit b ^ A, i.e. far away 
from the insulator-semimetal transition. In general, con- 
ductivities in all three directions are different. Note that, 
obviously, an externally applied field will add to the inter- 
nal exchange field b due to the ordered magnetic impurity 
moments, leading to a linear dependence of the conduc- 
tivity on the field. This would be supplemented by an 
orbital contribution, not discussed here. 

The optical conductivity of the line-node semimetal 
can be expected to behave as a function of frequency in 
the same way as the optical conductivity of graphene, i.e. 
to be roughly frequency-independent at low frequencies. 

As a final note, let us mention the expected behav- 
ior of the doped line-node semimetal. As in the case 
of the Weyl semimetal, discussed above, we will assume 
that the ionized dopants of density Ui (say donors) act 
as long-range Coulomb impurity scatterers for the doped 
carriers. Following the same line of reasoning as in the 
Weyl semimetal case, we then obtain: 



h o?b'^ 



(107) 



where we have ignored the Fermi velocity anisotropy for 
simplicity. The conductivity thus scales linearly with 
the carrier density, the same result as in graphene.'^Sl 
This similarity, however, is accidental in this case, as in 
graphene the linear scaling is obtained in a very different 
physical situtation: the Coulomb scatterers are charged 
impurities in the substrate, while the finite carrier den- 
sity is provided electrically by applying gate voltage. 



VI. DISCUSSION AND CONCLUSIONS 

In this paper we have considered two classes of topo- 
logical nodal semimetals, both of which occur in a multi- 
layer heterostructure, made of thin TI films, separated by 
ordinary-insulator spacers. Topologically-stable nodes, 
in which conduction and valence bands touch, occur in 
this system when TR symmetry is broken, by either mag- 
netic impurities or external magnetic field. Both point 
and line-node semimetals are characterized by protected 
surface states. These are especially robust in the case of 
the point-node, or Weyl, semimetal. The edge states in 
this case are chiral quantum Hall edge states, their chiral 
character making them robust even to hybridization of 
the bulk Dirac points.!^ 

The surface states of the line-node semimetal are "flat 
bands" , i.e. they are approximately dispersionless in a 



subset of the surface BZ, bounded by a projection of the 
bulk nodal line onto the surface plane. Since a flat band 
has a divergent density of states, nontrivial correlation 
effects, e.g. superconductivity or magnetism, may be ex- 
pected.i^ 

We have discussed transport properties of both types of 
nodal semimetals, as these can be expected to be impor- 
tant in the experimental characterization of these phases. 
We summarize those of the Weyl semimetal, including 
those not derived here, for convenience. If time-reversal 
symmetry is broken, it may exhibit an anomalous Hall 
conductivity. A general expression for this isP 



(108) 



where K is a wavevector, which can be expressed in terms 
of the Weyl points according to: 



K = Kn 



(109) 



Here qi = ±1 is the charge (in units of quantized U{\) 
Berry flux) of the Weyl point located at k = , and Kg 
is a reciprocal lattice vector (which could be zero). The 
former is a quantized anomalous Hall contribution due to 
completely filled bands. Note that although Eqs. (1081, 



( 109 1 generally describe a non-quantized anomalous Hall 



effect, it can be considered to be "semi-quantized" , in the 
sense that if K is measured experimentally, the universal 
quantized prefactor e^/h can be extracted. In our mul- 
tilayer case the Hall conductivity is non-zero and semi- 
quantized in this way even in the absence of an applied 
(orbital) field. In the proposed Weyl semimetal in the 
iridium pyrochlores,'^ it vanishes in zero applied field by 
cubic symmetry. However, an anomalously large Hall co- 
efficient may be induced according to Eqs. (1081, (1091, 



as the Weyl points shift in an applied magnetic field. 

In this paper, we showed that the bulk diagonal con- 
ductivity in the Weyl semimetal exhibits approximate 
quantum critical scaling due to Coulomb interactions. 
This implies that the zero temperature DC conductivity 

1/3 

is proportional to n,- , where rij is the density of charged 
donor impurities, and that the dependence of the conduc- 
tivity on temperature and frequency is approximately a 

1/3 

universal function of vpn-' /T and ui/T. We note that 
this conductivity scaling is a general property of any Weyl 
semimetal, including not only the superlattice structures, 
described here, but also bulk realizations, such as pro- 
posed for iridium pyrochlores.*^ A recent preprint draws 
similar conclusions in a model with rii = 

Finally, the diagonal conductivity also gets surface and 
interface contributions, due to edge states. In fact, even 
in a single crystal, Ising magnetic domains may form, and 
there can be chiral surface states, bound to such a domain 
wall. In the superlattice model of Sec. Ill A[ this is indeed 
the case for any domain wall, which is not normal to the 
z axis. If a sufficient density of such domain walls are 
present, they may give an appreciable contribution to 



16 



the diagonal conductivity. In practice, such contributions 
should be extracted by a careful study of hysteresis and 
by finding ways to align the magnetic order into a single 
domain. 

Assuming the line node to lie at the Fermi level (which 
is undoubtedly an approximation as discussed in depth in 
Sec. IVB 1 ), the resulting transport properties are some- 



what similar to graphene, except for the fact that these 
occur in a 3D material in our case. This is not unex- 
pected, as a line node in 3D is equivalent to a point node 
in 2D, since a point node can be thought of as a sec- 
tion of the line node by a plane in momentum space. As 
a consequence, a 3D line node has the same low-energy 
density of states as a 2D point node, i.e. linear in energy. 
This, in turn, leads to similarities in the transport prop- 
erties. In particular, the DC conductivity of the line- node 
semimetal is "universal" , in the sense of being indepen- 



dent of disorder. It is not, however, as universal as the 
conductivity of graphene, as it does depend on other ma- 
terial parameters, like Fermi velocities, the magnitude of 
spin splitting and the tunneling matrix elements, charac- 
terizing the TI and ordinary-insulator layers in the het- 
erostructure. 
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